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Two-Particle correlations based on the interference of identical particles has provided the chief 
means for determining the shape and lifetime of sources in relativistic heavy ion collisions. Here, 
Strong and Coulomb induced correlations are shown to provide equivalent information. 



In the collision of highly relativistic heavy ions, a 
highly excited region of matter is created which occupies 
thousands of fm 3 and explodes and dissolves on a time 
scale of 10 fm/c. Only the momenta of the collision debris 
are experimentally accessible. However, space-time infor- 
mation is crucial for understanding the reaction Q . Most 
importantly, a latent heat associated with a phase transi- 
tion would be accompanied by a significant enhancement 
to the duration of the emitting phase0, 0, 0, 0]. The 
most direct means for determining the spatial and tem- 
poral characteristics of the reaction is through the exper- 
imental measure of the two-particle correlation function 
which describes the ratio of the two-particle probabil- 
ity to the product of the single-particle probabilities for 
emitting species a and b with momenta p a and p\>- 

The correlation of two particles with outgoing mo- 
menta, p a and pb, is related to the probability of the 
source points being separated by a distance r, g(r), as 
measured in the two particle center-of-mass frame. 

C(q)-|d 3 r9(r)|^(q,r)| 2 . (1) 

Here, q = p„ = — P& is the relative momentum in the 
center-of mass frame, and g(r) describes the normalized 
probability that the source points of the two particles 
are separated by a distance r. If the particles do not 
interact with one another, the wave function, ^(q, r), 
would be e lq r , and the resulting correlation would be 
unity. Since the correlation function can be analyzed 
for any total momentum, p a + Pb, and for any choice of 
species a and b, a wealth of information regarding the 
space-time characteristics of the source is available. 

A measurement of C(q) represents three dimensions 
of information and can lead to a unique extraction of 
g(j). The temporal information contained in the four- 
dimensional source function S(r,t), which describes the 
probability of emitting particles separated by r and time 
t is folded into <?(r), 

g(r) = J dtS(T-vt,t), (2) 

where v is the average velocity of the emitted pair. If 
the characteristic time scale of S(r, t) is large, g(r) will 
be spread out along the direction of v relative to the other 



dimensions. In correlation analyses, it is standard to de- 
fine the direction of v as the "outwards" direction. The 
other two dimensions are referred to as "longitudinal" 
(along the beam axis) and "sidewards" (perpendicular to 
the beam and to the outwards direction) . Assuming that 
the spatial dimensions of S(r, t) are roughly equivalent in 
the outwards and sidewards direction, one can estimate 
the lifetime by the difference in the two dimensions of 
5(r), 

V 2 T 2 1=3 J? Qut — i?s idc , (3) 

where i? ou t and i? s id e are the spatial dimensions extracted 
by a three-dimensional analysis of C(q) with no consid- 
eration of lifetime. This estimate for r should work well 
for large lifetimes. 

Remarkably, analyses of two-pion interferometry at 
RHIC show equal outwards and sidewards sizes which 
suggests extremely sudden emission 0, 13 ■ Dynamical 
models that incorporate a phase transition with a latent 
heat strongly disagree with this observation 0, 0j . The 
two-pion analyses are based on the interference associ- 
ated with symmetrizing same-sign pions. As the exis- 
tence or non-existence of a latent heat with the QCD de- 
confinement transition represents a central issues of the 
RHIC program, it is imperative to explore alternative 
means for measuring the size and shape of the emission 
sources. 

For non-interacting identical particles the wave func- 
tion has a simple form, 

|V(q,r)| 2 - l±cos(2q-r). (4) 

By performing an inverse Fourier transform of C (q) — 1 in 
Eq. one can obtain g(r) . Even though the pions are 
charged, experimental analyses have tried to ignore the 
Coulomb interaction as much as possible, and in fact try 
to "correct" their data as to best eliminate the effects of 
Coulomb from the correlation function so that compari- 
son with simple forms for g(r) is easily accommodated. 
The aim of this paper is to demonstrate that Coulomb 
and strong interactions between the pair should not only 
be included in correlation analyses, but that they provide 
tremendous insight into both the size and shape of g(r). 

First, we consider two particles which interact only via 
the Coulomb interaction. Correlations for pK + from a 
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FIG. 1: pK + correlations are shown for a Gaussian source 
(-Riong = fisidc = 4 fm, i? ut = 8 fin). The classical approx- 
imation well explains Coulomb correlations at large relative 
momentum. The strong interaction only moderately affects 
the correlation function. 
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arc displayed in Fig. ^ as a function of Q inv = 2q. The 
source sizes are chosen i? ut = 8 fm and R± = 4 fm. 
The integration described in Eq. was performed with 
Monte Carlo methods. The Coulomb plane waves are 
solution to Schrodinger's equation for a wave where the 
outgoing wave has momentum q. 



-V 2 V(q, r) + q— V>(q, r) = q 2 ^, r), 
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where r\ is the Sommcrfcld parameter, 



V = 



Z a Z b ^ie 2 



EL EL 



9 



(E' a 



E' b ) 



(0) 



(7) 



The usual definition of the reduced mass /i, involving 
m a and nib, has been altered to achieve consistency with 
relativistic treatments for small e 2 . 

In order to understand the form of the squared wave 
function, it it insightful to compare to the classical ap- 
proximation where the squared wave function in Eq. Q 
is replaced by the ratio of the initial and final phase space. 



lV>(q,r)| 2 



d 3 q _ 1 + cos6> gr . - rj/(qr) 

<Pq + cos e qr ) 2 - 2(1 + cos O qr )ri/(qr) 

■Q(l + cos9 qr -2ri/(qr)), 

(8) 

where 8 qr is the angle between q and r. Integrating over 
cos 9 qr , one finds the angle-averaged correlation weight, 
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FIG. 2: pK + correlations are shown as a function of the angle 
of the relative momentum relative to the outwards direction 
for the Gaussian source (iiiong = -Rside = 4 fm, R ou t ~ 8 fm). 
The classical approximation becomes reasonable for large Qinv 
where the ratio of the suppression at cos 6 = to the suppres- 
sion at cos# = approaches (Rout / Reidc) 2 ■ 



In the non-relativistic limit this approaches unity as 1/q 2 . 
Convoluting this expression with the Gaussian source 
function gives a classical result for the correlation func- 
tion which is remarkably close to the quantum correlation 
function for large q as can be seen in Fig. ^ Thus, the 
tail of the correlation function provides a measure of the 
expectation, (1/r). 

Extracting the shape of the source requires measuring 
the correlation function as a function of the direction of 
q. Since the correlation approaches unity as 1/q 2 , we 
recommended plotting q 2 [C(q) — 1] rather than C(q) so 
that the main q dependence can be ignored, allowing the 
use of larger bins in q. 

As a function of the direction of q, correlations are 
shown in Fig. for Q inv = 30 and 150 MeV/c along- 
side the analogous results calculated with the classical 
form described in Eq. JSJ. The negative correlation at 
cos 9 = ±1 derives from the fact that the relative momen- 
tum of the two positive particles are deflected away from 
the direction defined by their relative position by the re- 
pulsive Coulomb force. Careful analysis of the classical 
expression, Eq. (JHJ, shows that for small 77, the nega- 
tive correlation is confined to q and r being anti-parallel. 
Since the chance of the particles being separated along 
the outwards axis is larger by a factor of R^/R 3 ^, and 
since the strength of the correlation in the two directions 
scales as R±/R ou t, the strength of the Coulomb induced 
correlation will be stronger along the i? Q ut direction by a 
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factor of i?out/^sidc at l ar S e 1- 

It should not be particularly difficult to obtain the re- 
quired statistics even though the correlation is of the or- 
der of 1%. Although the correlation falls as l/<2 2 nv , phase 
space rises as Qf nv . Unlike identical-particle correlations 
at small relative momentum, there are no issues with 
two-track resolution. However, one has to consider large- 
scale correlations, e.g., collective flow, jets and charge 
conservation. Correlation from charge conservation can 
be neglected by considering only same-sign pairs. Col- 
lective flow can be eliminated by carefully constructing 
the correlation denominator with pairs from events with 
the same reaction plane. Such competing correlations 
ultimately limit the range in Qmv that is useful for the 
the analysis. For instance, if the uncertainty of the com- 
peting correlations for pK + is of the order of 1%, the 
analysis should be restricted to Q; nv < 150 MeV/c so 
that Coulomb remains the dominant factor. 

The angular sensitivity of Coulomb-induced corre- 
lations have been investigated in intermediate-energy 
heavy-ion collisions. These studies involved light frag- 
ments, e.g. Carbon-Carbon, which can be treated with 
classical trajectories The Coulomb mean field from 
the remainder of the source was found to distort the the 
shape information. However, the residual Coulomb inter- 
action is expected to play a much smaller role at RHIC 
where the hadrons move much faster and spend less time 
interacting with the mean field. Furthermore, analyses 
can be performed with both positive and negative pairs, 
e.g., both pK + and pK~ . Averaging the two correlation 
functions should largely cancel the effects of the resid- 
ual interaction. At high energy, the angular dependence 
of Coulomb correlations of non-identical particles have 
been used to determine the degree to which one species 
is displaced relative to another. These studies involved 
comparing correlations for 9 < 7r/2 with correlations with 

0>tt/2 Eim. 

Incorporating the strong interaction into the wave 
function in Eq. is straight-forward if the the strong in- 
teraction can be expressed as a non-relativistic potential. 
One must simply solve the Schrodinger equation. How- 
ever, many interactions can not be expressed in terms of 
such potentials. For instance, the interaction of a n + and 
proton through the delta resonance involves a quantum 
rearrangement of the participating quarks. Thus, a wave 
function is not a meaningful quantity for separations be- 
low e ~ 1 fm, but it is certainly a well defined object for 
r > e, and can be expressed in terms of a Coulomb wave 
with the incoming partial waves are modified by phase 
shifts. 



Here, Fi and Gi are the regular and irregular partial 
Coulomb waves and the second term describes the dis- 
tortion of the incoming partial wave, Fg — iGe- 

For the pir + and pK + examples discussed in this study, 
the plane wave also has a spin index, m s = ±1/2, which 
is not conserved for partial waves with I > 0. The states 
of a given m s must be decomposed in terms of eigenstates 
of total angular momentum which are phase shifted by 
eigen-phases, 5j^. After applying some angular momen- 
tum algebra, Eq. ifTUI) can be modified to include flipping 
the spin. For I = 1, 
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Thus, the wave function for r > e can be evaluated if 
given the phase shifts. 

For r < e, one must use an effective form for r)| 2 . 
Since e is much smaller than any characteristic dimension 
of the source, only the integral of ip 2 matters in the region 
less than e, and one can safely choose, 



|^(q,r< e )| 2 = |^o(q,r)| 2 + iy( e ,g), 



(12) 



where ipo is the Coulomb wave and W is independent of 
r. The change in the density of states can be ex pres sed 
both in terms of phase shifts and wave functions |l2j , 
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Thus, W can be expressed in terms of derivatives of the 
phase shifts and integrals of the type, 
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Since <pe is a solution to the Schrodinger equation, one 
can rewrite if, assuming that <fi and </>' are solutions with 
eigenvalues q and q' ~ q, 



V(q, r) = </>o(q, r) + £ y/4n(2£+l) 
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FIG. 3: pix + angular correlations for a Gaussian source 
(-Rout = Rsidc = 4 fm, 7?i ong = 8 fm) with Q inv = 450 MeV/c. 
At this Qinv the correlation is dominated by the A ++ reso- 
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±1 can be explained by 



Transforming the derivatives to the variables rj and 
x = qr facilitates use of the recursion relations for the 
Coulomb wave functions [l3| . 
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The upper limit, x = oo, need not be evaluated since 
it would be canceled by an equal and opposite contribu- 
tion from the non-phaseshifted term in Eq. (|13|l . For the 
case with no Coulomb interactions, one can derive a sim- 
pler form involving spherical Bessel functions rather than 
Coulomb wave functions. For the no-Coulomb example, 
one can show that W — as q — ► 0. 

Using experimentally tabulated phase shifts, and the 
arbitrary choice of e = 1 fm, wave functions were numer- 
ically generated and convoluted with the source function 
through Eq. (JIJ to generate correlation functions. Corre- 
lations for pK + were only slightly affected by the strong 
interaction as can be seen in Fig.s ^ and [21 For this 



example, there are no resonant channels and the corre- 
lation is dominated by Coulomb. Figure shows results 
for pir + correlations which are dominated by the A ++ 
for Qi nv near the resonant momentum, Qi nv = 450 MeV. 
A strong angular correlation is caused by shadowing in 
the forward direction. Like Coulomb-induced correlation, 
there is a dip when the relative momentum is aligned with 
the long-axis of the source. 

Strong and Coulomb induced correlations had been 
previously studied for their ability to to unfold the angle- 
averaged source function, q(r) , for both high-energy and 
low-energy collisions [141 Il5l fl6| . Our findings show 
that such correlations also have tremendous potential to 
discern shape characteristics and thus provide an esti- 
mate of source lifetimes. Determining shape and life- 
time characteristics had been previously confined to anal- 
yses of identical-particle correlations. As strong and 
Coulomb correlations are of a manifestly different charac- 
ter than correlations from identical-particle interference, 
the analyses described here represent a truly indepen- 
dent strategy for determining space-time characteristics 
of hadronic sources. 
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